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ABSTRACT 



In this paper we present the tools to optimally extract the Lensing-Integrated Sachs Wolfe (L-ISW) bispectrum signal from future 
CMB data. We implement two different methods to simulate the non-Gaussian CMB maps with the L-ISW signal: a non-perturbative 
method based on the FLINTS lensing code and the separable mode expansion method. We implement the Komatsu, Spergel and 
Wandelt (KSW) optimal estimator analysis for the Lensing-ISW bispectrum and we test it on the non-Gaussian simulations in the 
case of a realistic CMB experimental settings with an inhomogeneous sky coverage. We show that the estimator approaches the 
Cramer-Rao bound and that Wiener filtering the L-ISW simulations gives a slight improvement on the estimate of /j^^" "' of < 10%. 
For a realistic CMB experimental setting accounting for anisotropic noise and masked sky, we show that the linear term of the 
estimator is highly coiTelated to the cubic term and it is necessary to recover the signal and the optimal error bars. We also show that 
the L-ISW bispectrum, if not correctly accounted for, yields an underestimation of the /^"°' error bars of - 4%. A joint analysis of the 
non-Gaussian shapes and/or L-ISW template subtraction is needed in order to recover unbiased results of the primordial non-Gaussian 
signal from ongoing and future CMB experiments. 

Key words. The Cosmic Microwave Background, non-Gaussianity, Lensing, ISW, cosmology 



1. Introduction 

One of the most relevant mechanisms that can generate non- 
Gaussianity from secondary Cosmic Microwave Background 
(CMB) anisotropics is the coupling between weak lensing and 
the Integrated Sachs Wo lfe (ISW) (ISachs & NVolfd 119671) and 
the Rees Sciama (RS) (lRees&Sciamalll968h . This correla- 
tion gives in fact the leading contribution to the CMB sec- 
ondary b ispectrum w ith a blackbody frequency dependence 
(iGoldberg & Scergel 19991 IVerde & Sp ergel 2002; Giovi et aL 
120051) . Weak lensing of the CMB is caused by gradients in 
the matter gravitational potential that distorts the CMB pho- 
ton geodesies. The ISW and the RS effects, on the other hand, 
are related to the time variation of the gravitational potential 
wells. The relevant mechanism is given by the late ISW, ow- 
ing to the action of Dark Energy which causes the decay of 
the gravitational potential wells as the Universe expands. Both 
the lensing and the ISW effect are then related to the matter 
gravitational potential and thus are correlated phenomena. This 
gives rise to a non-vanishing three-point correlation function 
or, analogously, a non-vanishing bispectrum, its Fourier coun- 
terpart. The RS (also referred as the non-linear ISW) arises 
when the growth of structure in the evolving universe becomes 
non-linear. Being a second order effect, the RS gives a smaller 
contribution to the signal with respect to the ISW. The CMB 
bispectrum arising from the cross correlation between lensing 
and ISW/RS (from now on referred to as L-ISW) is expected 
to have an high signal-to-noise ratio from ongoing and future 
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CMB experiments so that it will be de tectable in the near fu- 
ture with an high statistical significance (Verde & SDerge lll2002l : 
Giovi et al. 2005; Mangilli & Verde 2009; Lewis et al. 201^ A 
detection would open the possibility to exploit the cosmolog- 
ical information related to the late time evolution encoded in 
the L-ISW signal. It is useful to stress that a significant de- 
tection of the L-ISW signal from ongoing CMB experiments 
like Planck would be a powerful probe of Dark Energy from 
CMB alone and it would be a complementary probe of the late 
time Universe with respect to the large scale structure and the 
the CMB power spectrum an alysis. Moreover. iMangilli & Verdel 
( l2009h ; lHanson etalj ( |2009|) showed that the L-ISW bispectrum 
can be a serious contaminant problem for the estimation of the 
primary local non-Gaussiani ty from future da ta. Ongoing CMB 
experiment such a s Planck ( Ade et all 201 ll) and future exper- 
iments like COrE (iBouchet et al. 201 1 ) will then require a de- 
tailed reconstruction of the L-ISW bispectrum either to be able 
to correctly remove the L-ISW contribution when estimating the 
local primary non-Gaussian parameter fyi, or to exploit the cos- 
mological information encoded in the signal; therefore it be- 
comes extremely important to know how to model and simulate 
it. 

In this paper we present the formalism and the numerical 
implementation /) to generate simulated CMB maps containing 
the L-ISW signal and //) to build and test the optimal estimator 
for the L-ISW bispectrum, accounting for both the cubic and the 
linear parts. The linear pait for this specific kind of signal, has 
been here calculated and tested for the first time. As regarding 
the CMB non-Gaussian simulations, we implemented and tested 
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the L-ISW si gnal with two methods : the separable mode expan- 
sion method ("Fergusso n et al.lQoiOl: ISmith & Zaldarria"gall201 ih 
and the non-perturbative approach described in Sec. 12.1] 

It is important to have an optimal estimator for the L-ISW 
bispectrum in order to extract the signal optimally from future 
data and to disentangle it from other kinds of non-Gaussianities, 
i.e. the local primar y bispect r um, w ith whi ch it is degener - 
ated. Here, following lKomatsul <2010l) and iMunshi et alj (1201 Ih. 
we im plemented the KS W bispectrum estimator ( Komatsu et alj 
l2005h for the L-ISW signal of a full sky, cosmic variance lim- 
ited CMB experiment and in the case of a more realistic in- 
strumental setting, similar to that of a space-based expeiiment. 
Furthermore, for this realistic case, we investigate the statistical 
detection significance and the impact that the L-ISW bispectrum 
has on the estimation and on the variance of the primary local 
non-Gaussian parameter f^i- 

The outline of this paper is as follows. In section|2]we present 
the methods to simulate the non-Gaussian CMB maps contain- 
ing the L-ISW bispectrum signal by the use of both the separa- 
ble mode expansion method and the non-perturbative covariance 
method. Section[3]provides the basics to build and implement the 
optimal estimator for the L-ISW signal, including its linear part. 
It also includes a discussion regarding the implementation of the 
Wiener filtered simulations algorithm. In section |4] we present 
the relevant tests and results. In section |5] we quantify the sta- 
tistical detection significance of the L-ISW bispectrum and the 
impact on the error of primary non-Gaussianity f^i due to the 
presence of the ISW signal. Finally, in section|6] we discuss the 
results and we summarize the conclusions. Details on the simu- 
lations built with the covariance method and on the L-ISW cross 
coiTelation coefficients are given in the appendix. 



2. Simulated non-Gaussian CIVIB maps 

In this section, we present the formalism to create simulated 
CMB maps for the L-ISW bispectrum. We use two differ- 
ent methods: a no perturbative approach, here named the 'co- 
variance method' , and the se parable modes expan sion method 
(( iFergusson etal.ll20I0h and dSmith & Zaldarriagall20TTl) ). The 
latter gives an efficient and easy to handle way to generate L- 
ISW maps, while the former method provides better insights on 
the physics related to the L-ISW bispectrum. In this case, in fact, 
the L-ISW signal is generated starting from the covariance ma- 
trix representing the expected correlation between the lensing 
and the ISW/RS effects. 



2.1. Covariance method 

The L-ISW correlation is defined by the covariance matrix: 



ISW 



cr c 

f~.T<t> J 



(1) 



(the gravitational potential projection along the line of sight) is 
defined by: 



(n) = -2 dr 0(r, nr). 

Jo nz) r(zis) 



(3) 



The term in the numerator, , is the power spectrum of the 
cross correlation between the lensing and the ISW/RS effect, see 
appendixlBlfor details. 

After a Cholesky decomposition of the L-ISW correlation 
matrix Cl-/sw, the two new variables ff„, and zcm are then de- 
fined by 



Zlm 




+ y{,. 



(4) 
(5) 



where X(„, and yf„, are two independent random gaussian fields. 



By definition, the new fields are such that: {t~} = 
C^J and they have the non-zero cross correlation {zt) - C 



As described in appendix |A] the map that contains the de- 
sired L-ISW bispectrum is then given by the coefficients 

= ze„, + flL - 4n = ze.n + Afl^„,, (6) 

where a^^^^ and a^^^^ are, respectively, the unlensed primary and the 
lensing angular coefficients and Aajr^^ - a^^^ - a^^^ corresponds to 
the lensing expansion terms only. Note that by construction y(,„ 

has the same phases as a^^^ = yc„, yjcf^ and 0^^^ = X{,„ -^Jc^ 

the same as Xf,„, which is necessary for building a map with the 
wanted bispectrum signal. 

Figure [1] shows in black the temperature CMB power spec- 
trum of one simulated L-ISW map, C^^'^^, built from Eq. |6] 
The non-Gaussian contribution, in blue in the figure, is always 



subdominant and the C 



L-ISW 



are consistent with the theoreti- 



cal input (Cj^)ih (red line) obtained with CAMB dLewis et al.l 



from map 
input from CAMB 
C,™ =(C,""-C^) /rom map 




1500 



and the cross correlation coefficient: 



C 



Here. 



(2) 



are respectively the CMB primary temperature power spectrum 
and the lensing power spectrum, where the lensing potential 



Fig. 1. The L-ISW power spectrum from the covariance method 
simulation. The plot shows that the temperature power spectrum 
of the L-ISW simulations generated with the method described 
in Sec. 12. H is compatible with the input theoretical power spec- 
trum from CAMB and that the non-Gaussian contribution is al- 
ways subdominant. The temperature power spectrum from one 
simulated L-ISW realization is shown in black, the red line refers 
to the theoretical input from CAMB while the blue refers to the 
non-Gaussian L-ISW contribution from the same realization. 
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and simulations of the the CMB Lensing-ISW bispectrum 
Here, 

Pin) = ^a{„,Y(,„(n), 

Cm 

Q{n) = Y^C]\C-'a)c,nYUn), 

Cm 



1500 



Fig. 2. The L-ISW power spectrum from the separable mode ex- 
pansion method simulation. The temperature power spectrum of 
the L-ISW simulation generated with the method described in 
Sec. l2.2l is compatible with the input theoretical power spectrum 
from CAMB and the non-Gaussian contribution is always sub- 
dominant. The temperature power spectrum from one simulated 
L-ISW realization is shown in black, the red line refers to the 
theoretical input from CAMB while the blue refers to the non- 
Gaussian L-ISW contribution from the same realization. 



120001) B As throughout the paper, the reference cosmological 
m odel used is the ACD M model with parameter values defined 
in (iKomatsuet alJl201lh . 

2.2. Separable modes expansion method 

Follow ing iFergusson et al.l ( 1201 Ol) and ISmith & Zaldarria"gal 
(1201 lb . the non-Gaussian part of the CMB angular coefficients 
can be defined starting from a given reduced bispectrum. In the 
case of the L-ISW signal, the method can be used because this 
kind of signal is separable, so 



ycitm iL-ISW 



3? yf(n)aG Y';'(n) 



(7) 



From the expression of the L-ISW reduced bispectrum in Eq. 
( fT3] l and by factorizing the I dependence, the explicit form of 
the non-Gaussian contribution to the Of,,, from the L-ISW cross 
correlation is given by; 



Vatm \l-isw 



d^nY'"(n) 



^(^+l)e(n)£(n) 



+ Q([52£](n)Q(n) - [s'QmEin)) 

- ([5V](n)e(n) + [6'-QmP{n)) 

- {{{ + l)Q{n)P{n) 

+ qc[[6'EmP{n) - [6^P](n)E(n)) 

+ {({ + l)q(Pin)Ein) ■ 



(8) 



http://camb.info 



(9) 



Cm 



The maps with a prefix are given by, e.g., S^P - - + 
l)a(inY(in{n); they correspond to the maps of Eq. [TOl multiplied 
by the -{{{ + I) factor. The final solution containing the L-ISW 
signal is then: 



(10) 



where a^^^ is the Gaussian part. 

In Fig.|2] we show the CMB temperature power spectra from 
the Gaussian and the non-Gaussian map, as defined in Eq. dS). 
The non-Gaussian contribution is always subdominant as ex- 
pected. 

3. The Optimal KSW estimator for the 
lensing-ISW/RS bispectrum 

In this s ection we present the formalism related to the KSW es- 
timator (iKomatsu et alJl2005l) for the Lensing-Integrated-Sachs 
Wolfe bispectrum signal. 

3.1. Definition 

The ac„, probability distribution function (PDF) in the limit of 
weak non-Gaussianity (i.e. tmncated at the bispectmm level) is 
given by (,Babich,2005.: .Tavlor & Watts 2001; Komatsu. 20 10.) : 



P(a) 



1 



■ exp 



1 all Ijifij 



1 



V V al(C-\, 



-2 Zulu""" 

Im I'm' 



(27r)^i,™/2|C|i/2 

^l\m\ ^hmj^hm^ 

)[(C ^a)i^m^(C 'fl);2m,(C '«)/„„, 



(C '«)/„„,]}■ 



(11) 



where {a(^„,^a[^j„^_a(^,„^} is the angular bispectrum. Here, we are 
interested in the L-ISW case, for which the angular bispectrum. 



parametrized by the amplitude parameter f^j^^^ , is 



where 



,m,m,Bi, rL-ISW , L-ISW 



(12) 



+ 1) - hih + 1) + hih + 1) , 
7 C{Cf; + (5p) 



,(13) 



is the reduced bispectrum and C^*^ = {(t>*[...a^7j ^'^) are the L- 
ISW cross-correlation coefficients. According to ( iKomatsu et al.l 
l2005h . for small departure from Gaussianity, the optimal estima- 
tor for theL-ISW amplitude parameter is given by: 



JNL - )'J L-ISW, 

where (F ') is the inverse of the L-ISW Fisher matrix 



nL-ISW nL-ISW 

tit2[i hhh 



(14) 



(15) 
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Fig. 3. The Large scale contribution to the non-Gaussian L-ISW 
signal. Upper panel: the map g(n) = Z(m C^/{C'^a)(„Y(,n{n) 
contains the L-ISW coefficients C^^ and enters the L-ISW esti- 

mator Eq. [18] The ^-filter acts as a filter which suppresses 
the small scales (lower panel). 



Fig. 4. Tlte intermediate scale contribution to the non-Gaussian 
L-ISW signal. Same as figure|3]but for the map ^-Y.ciii + 
l)C^'^(C^^a)(mYcm{^) (upper panel) and its corresponding filter 



-H 1)^. The factor 



+ 1) dominates at high { defining 
more small scale features with respect to the previous Q map. 



In the case of a realistic CMB experimental setting, the noise, 
A';, and the beam window function, W(, are accounted for so that 
Ce =N(-^- Cfw]. In this case, the bispectrum is also convolved 
with the beam transfer function wc, B^^^'f^^ oc b'^^^'f^^ w^we^we-^. 
Given a mask M{p), the observed sky fraction fsi^y is defined as: 



Iskv 



N, 



(16) 



where Npix - 12N^ is the number of pixels in the map, A^, is the 
map resolution and the sum J]p is done over the pixels. 

Assuming that the only relevant non-Gaussian contribution 
is coming from the L-ISW term, which is the case if the local 
primordial non-Gaussianity is small and foregrounds and point 
sources have been correctly removed and masked, Sl-isw is 
given by the data as 



all Im 



(C 'a)f,,„,(C 'fl)f,,„,(C 'a)f3„„ 



3(C ')c,mi,f2m2(^ 



By factorizing the £i dependence, this becomes 



(17) 



> L-ISW 



1 J d^n{p(n 



H- S 



lin 



(18) 



where the maps P(n), ^(n), Q(n) etc. are the same as defined 
in Eqs. \W\ and, in the case of a realistic experiment, they are 
convolved with the experimental window function W( so that, 
for example, P(n) s Y.Cm wcac,„Y(,„{h). 

In Eq. [18] the first two lines refer to the cubic part of the 
estimator, while S^^'^^ is the linear part which corrects for 
anisotropics and must be included in the case rotational invari- 
ance is not preserved. Details on the analytic expression of the L- 
ISW linear term and on its numerical implementation are given, 
respectively, in the next subsection l3.2l and in Sec. [4] 
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Fig. 5. The plot shows the f^^^'^^ values for 100 simulated non- 
Gaussian maps obtained with the covariance method of Sec. 
12.11 The lensing part has b een computed with the FLINTS code 
dLavaux & Wandeltll20ldl) . The straight line refers to the aver- 
aged f^2'^^ from these simulations, while the dashed line to the 
averaged 1-cr error Here {„uix - 1000. 



3.2. The linear term 

The linear term of the estimator is given by 



aL-lSW 1 r ^2,- 'V hL-IS\ 

all €m 



(19) 



{c-')t,m, (c-'fl),„„3 y;;- m';; (fl)^"' (ft) 

ng the explicit 1 
f-dependence one obtains 



By using the explicit form of b'^J^i'^ and by factorizing the 



aL-ISW 



= _i J cfh\Q{h)Yp{n)5'-E{h))Mc 

- 6^Q(h){P(h)E(mMC 

- E(h) 

+ 6^E(h){P{fi)Q(h))Mc 

+ P{h) 



- {E{h)6^P{h))Mc 
(20) 

{Q{h)5^P(h))Mc - {Pin)S~Qin))MC 

- 6^P(h){E(h)Q(h)}Mc 



{Q(h)6^E(h)}MC - {E{h)6'Q{h))Mc 



where {)mc indicates the Monte Carlo (MC) averages and the 
different maps are defined in Eq. [10] and they are convolved 
with the experimental window function Wf, so that P(n) = 
'Zem weatmYeJh), etc. 

3.3. Wiener filtered maps 

The optimal bispectrum estimator as described in Eqs. ([Tt] [TsT i 
involves products of inverse variance filtered maps, C 'a - 
{S + N)~^a, where S and are the signal and the noise co- 
variance matrix, respectively. A brute force calculation of such 
an expression is impractical for modern high-resolution exper- 
iments as it involves the inversion of two matrices that are too 
large to be stored and processed as dense systems. In case the 
noise covariance can be described in terms of a simple power 
spectrum in spherical harmonic space, the calculation simplifies 




Fig. 6. Same as Fig.|5]but for 100 simulations built with the sep- 
arable mode expansion method (Eq. (jSjl). 



significantly. However, this approach is no longer exact for ex- 
periments with anisotropic noise distribution or reduced sky cov- 
erage, leading to an increase in the error bars of the estimates. 

Here, we use Wiener filtering as a basis for the exact eval- 
uat ion of terms involving C 'a. We apply the iterative scheme 
of jElsner & Wandemi201 3l) to calculate the Wiener filter a^^ = 
S{S + Ny^a, the maximum a posteriori solution in case signal 
and noise are Gaussian random fields. After a^^ has been suc- 
cessfully computed, we finally obtain the inverse variance fil- 
tered map by normalizing the spherical harmonic coefficients of 
the Wiener filter solution by the CMB power spectrum multi- 



plied with the beam window function, C aim - (Cfbj) 



r/iA2\-l^WF 



4. Results 

In this section we present the results regarding the numerical 
implementation of the optimal estimator and of the methods pre- 
sented in sec. l2.1l and l2.2l to build the CMB maps containing the 
L-ISW bispectrum. In particular, we processed the simulated L- 
ISW maps through the estimator pipeline to get the amplitude 
parameter f^^'^'^ of Eq.[T4] We consider two main settings: 

- a full sky cosmic variance limited CMB experiment up to a 
maximum multipole {„,ax - 1000 and 

- a more realistic experimental setting which consists of a 
one channel CMB experiment with a Gaussian beam with 
a FWHM 6i, - 7', a galactic mask leaving ^ 80% of the sky 
and anisotropic uncorrelated noise. These settings are visu- 
alized in Figs. lC.3 l lC.2l and details are given in Sec.lCl 

All runs have been performed at full resolution Nside - 2048 
(which corresponds to a map pixel number of 5.033 ■ 10^). The 
map s in Eqs. (fTOli are calculated by using the Healpix pack- 
age (iGorski et al.ll2005l) . The theoretical power spectrum of the 
temperature-only primary CMB coefficients Cc has been gener- 
ated with the CAMB code for a fiducial ACDM cosmological 
model with parameters correspon ding to WMAP7 cosmological 
parameters (iKomatsu et al.ll20Tl1) . For illustrative purpose, the 
plots of the maps and of the correspondent ^-filters containing 
the L-ISW cross correlation coefficients C^'^ are shown in Figs. 
|3]and|4] 

We built a set of 100 CMB simulations for each of the two 
methods described in sections ISTSlandlTTIfor a cosmic variance 



5 



Mangilli, Wandelt, Eisner & Liguori: Optimal estimator and simulations of the the CMB Lensing-ISW bispectrum 



31. 



Z - 




20 40 60 80 100 



Fig. 7. The same as Fig|6]but for a more realistic CMB exper- 
iment with a 7' FWHM Gaussian beam, anisotropic noise and 
20% galactic mask. Here j. = 1500. The dashed lines are the 
1-cr averaged error bars from simulations while the dotted lines 
are the expected Fisher errors. 
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Fig. 8. The linear term of the estimator reduces the error bars 
in the case of anisotropic data. The plot shows the linear (solid 
black line) and the cubic (dot-dashed black line) contributions 
to the total (red line) /^^"'^ = (/^^^^ "'),„„/. + (/^^'^'^ )/,„«„• in 
the case of a CMB experiment with anisotropic noise and 22% 
galactic mask. 



limited CMB experiment with full sky co verage. For the co- 
varian ce method, we used the FLINTS code (iLavaux & Wandelg 
1201 oh to generate the lensing coefficients a'^^ and the lensing po- 
tential coefficients (pcm needed to build the non-Gaussian a^^'^^ 
as described in Sec. 12.11 In both cases, we analyzed the L- 
ISW simulated CMB maps with the L-ISW estimator up to 
fmax = 1000. According to the definition of f^l^^^ , the expected 
value is 1 with l-cr error predicted from theory for {„,ax - 1000 
of ^ 0.64. In the case of the separable expansions mode method, 
the simulations give a mean fj:^^'^^ - 1.1 with averaged 1- 
cr error ^ 0.69. With the simulations built with the covariance 
method, we obtain a mean fj:j^'^^ - 1-21 with averaged l-cr 
error of 0.67. The results are summarized in Fig. 0and Fig. 
respectively. These estimates are compatible with the theoreti- 
cal predictions. The error bars are slightly suboptimal because 
of numerical noise and the fact that we are assuming a diagonal 
covariance matrix so that (C^^a)i„, = af„JCf. 

In order to test the estimator on a more realistic case, we 
built a set of 100 simulations with the separable mode expansion 
method considering a realistic experimental setting. This con- 
sists of a CMB one channel experiment with a Gaussian beam 
FWHM 9h = 7', a galactic mask with /j^y = 0.78 and anisotropic 
noise, as previously described. In this case, we run the estimator 
up to Umax - 1500. The expected theoretical l-cr error on fj:^^'^^ 
for this experimental setting and up to C,„ax — 1500 is ^ 0.49. 
This estimate accounts for a ^ 10% percent increase in the error 
bar due to the fact that the lensing is intrinsically non-Gaussian 
and it gives an ex tra contribution to the variance, as shown in 
(iLewis et al.ll201lh . We get a mean /^'^'^ = 1.09 with averaged 
l-cr error ^ 0.55. In this case we computed both the cubic and 
linear part of the estimator In particular, the linear term has been 
tested with a set of 100 Monte Carlo (MC) averages generated 
for each map product in equation (Eq. ISTT i. In the presence of 
anisotropic noise and a sky cut, the linear part of the estimator is 
necessary to recover the expected estimation of f^^'^ ^ and error 
bars. The linear contribution to f^^^'^ ^ is strongly anti-correlated 
with the cubic part. This behavior is summarized in figure Fig. [8] 
In the plot are shown the linear and the cubic contributions to the 
total amplitude /^^"'^ = {f^-^^'^'^Ubic + iflii,"'^hnear- We also 



checked that with 100 MC averages the linear term converges 
and it is stable: for this specific experimental setting we find that 
the results do not improve when increasing the MC averages to 
200. 

Finally, in order to test optimality, we Wiener filtered the 
100 L-ISW simulations and we processed them through the L- 
ISW estimator pipeline. T he maps has been produced following 
(lElsner& Wandellll2013h . as described in section Sec. 13.31 We 
use as inputs the same experimental settings as described pre- 
viously. The linear term has been computed with 100 Wiener 
filtered MC simulations. We found that the improvement with 
respect to the non Wiener filtered simulations is small (< 10%) 
in the case of our particular settings. However, this does not ex- 
clude that the Wiener filtering may have a more noticeable im- 
pact for a more realistic experimental setting and noise covari- 
ance. 



5. /nl error estimation 

This section summarizes the results regarding the impact of the 
L-ISW signal on the error estimation of fyi from the local type 
non-Gaussianity. If the only contribution to fyi were from the 
primary local type non-Gaussianity the error on this parameter 
would be simply given by 



(21) 



i.e. the inverse of the Fisher matrix of the local type non-gaussian 
contribution 



(22) 



where fsky refers to the observed sky fraction. The noise, Nc, and 
the beam, b[, can be accounted for so that Cc - Nt + Cfb^. In 
this case the bispectrum is also convolved with the beam transfer 
function fef: B^^^.^^^ oc b^^^^^^be^bt^bt,. 
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However, the L-ISW can be a serious contaminant of the lo - 
cal primary signal (iMangilh & Verdel2009UHanson et al.ll2009h . 
so that it is important to quantify the effect on the expected fyi 
error as well. If the L-ISW signal is present, the error matrix will 
be given by a non-diagonal Fisher matrix of the form 



where 



- fsky ^ 



qL-ISW gP 



(23) 



(24) 



is the cross correlation term and F^^'^ ^ is the the Fisher term of 
the L-ISW signal of Eq. [15] The expected eiTor on the local ff4L 
will be then: 



(25) 



i.e. the inverse of the full Fisher matrix containing the cross cor- 
relation between the primary local non-Gaussianity and the L- 
ISW signal. The difference between the error estimation on f^i- 
primary with and without the L-ISW contribution is: 



A P P 



(26) 



To quantify the level of correlation between the two signals one 
can define the correlation coefficient as 



(27) 



We find that the effect of the local non-Gaussianity on the L-ISW 
is negligible. Therefore, the 1 - cr eiTor of the L-ISW amplitude 



parameter Z^;'"^ is given by cr^ 



For a realistic CMB experiment as described in Sec. |4]and 
appendix |C] we find that the correlation between the two sig- 
nals is r = 0.20 at l,nax = 1500 and r = 0.27 at (,„ax = 2000 
and that the expected detection significance of the L-ISW sig- 
nal {Ijcr^^'^^) is, respectively, at ^ 2 and ^ 3 cr. The effect on 
the ff^i local error due to the contamination is in the range be- 
tween ^ 3% to ^ 5 % for e,„ax from 1000 to 2000, depending on 
how much the two signal are correlated. For C„,ax - 2000, if the 
L-ISW signal is not accounted correctly, the fyi error bars are 
overestimated by ^ 4%. 



6. Discussion and conclusions 

We have presented the formalism and the numerical implemen- 
tation to build the optimal KSW estimator for the Lensing-ISW 
bispectrum. Moreover, we have tested the estimator on simu- 
lated CMB maps containing the Lensing-Integrated Sachs Wolfe 
(L-ISW) non-Gaussian signal and on the Wiener filtered simu- 
lations in order to test optimality. As regarding simulations, we 
have implemented and tested two methods: a non-perturbative 
approach to simulate CMB sky maps with the L-ISW signal 
which is based on the FLINTS lensing code dLava ux & Wandelll 
1201 Ol) and the perturbative separable mode expansion method 
calculated for this specific signal. We provide the analytical ex- 
pression and the numerical implementation of the linear term of 
the estimator for this specific kind of bispectrum. For a realistic 
CMB experimental setting accounting for anisotropic noise and 
masked sky, the linear term gives a relevant contribution which 
is highly anti-correlated with the cubic part and it is necessary 
to recover the signal and optimal error bars. In order to achieve 



optimality, we also tested the estimator on the Wiener filtered 
L-ISW simulated CMB maps. In this case we recovered the sig- 
nal with error bars which saturate the theoretical Cramer-Rao 
bound, with a small improvement of < 10% with respect to the 
non Wiener filtered simulations. Finally, we estimate that, if not 
correctly accounted for, the L-ISW effect has also an impact on 
/^™' error bars leading to a bias and an overestimation of ^ 4%, 
in agreement with (Lewis et al. 20X1)- Thus a joint analysis of 
non-Gaussian shapes and/or L-ISW template subtraction will be 
needed in order to recover unbiased minimum variance results 
of the local type primordial non-Gaussian signal. 

It is important to note that the KSW bispectmm approach 
to the estimation of the L-ISW is complementary to the lens- 
ing reconstruction estimator of (iLewis et al. I I2OIII) . In princi- 
ple, the KSW estimator can offer advantages with respect to 
other methods because the bispectrum has a unique shape and 
it has been shown to b e robust to foreground contamination 
jYadav & Wandellll2010h so it can be measured by using a larger 
sky fraction. In addition, inclusion of the L-ISW in the frame- 
work of bispectrum analysis gives an unified approach to testing 
for primordial non-Gaussianity. The tools presented in this paper 
enable the optimal analysis of this important signal from future 
CMB data. 

Acknowledgements. This work was supported in part by NSF grants AST 07- 
08849 and AST 09-08902, and by NASA/JPL subcontract 1413479; and through 
Ben Wandelt's ANR Chaire d'excellence ANR-lO-CEXC-004-01. AM acknowl- 
edges Guilhem Lavaux for the FLINTS lensing simulations, Licia Verde for use- 
ful comments and discussion and the University of Illinois for the use of the 
curvaton computers. 



References 

Ade, P. et al. 20 11, Astron.Astrophys., 536, 16464 
Babich, D. 2005. Phys.Rev., D72, 043003 
Bouchet, F. et al. 2011 

Eisner F & Wandelt, B. D. 2013, Astron.Astrophys., 549, Al 11 
Fergusson, J., Liguori, M., & Shellard, E. 2010, Phys.Rev., D82, 023502 
Giovi, F, Baccigalupi, C, & Perrotta, F 2005, Phys.Rev., D71, 103009 
Goldberg, D. M. & Spergel, D. N. 1999, Phys.Rev., D59, 103002 
Gorski, K., Hivon, E., Banday, A., et al. 2005, Astrophys.J., 622, 759 
Hanson, D., Smith, K. M., Challinor A., & Liguori, M. 2009, Phys.Rev., D80, 
083004 

Junk, V. & Komatsu, E. 2012, Phys.Rev., D85, 123524 
Komatsu, E. 2010, Class.Quant.Grav., 27, 124010 

Komatsu, E., Spergel, D. N., & Wandelt, B. D. 2005, Astrophys.J., 634, 14 
Komatsu, E. et al. 201 1, Astrophys.J.Suppl., 192, 18 
Lavaux, G. & Wandeh, B. D. 2010, Astrophys.J.Suppl., 191, 32 
Lewis, A. 2012, JCAR 1206, 023 

Lewis, A., Challinor A., & Hanson, D. 201 1, JCAR 1 103, 018 
Lewis, A., Challinor A., & Lasenby, A. 2000, Astrophys.J., 538, 473 
Mangilli, A. & Verde, L. 2009, Phys.Rev., D80, 123007 

Munshi, D., Valageas, P., Cooray, A., & Heavens, A. 2011, 

Mon.Not.Roy.Astron.Soc, 414, 3173 
Rees, M. J. & Sciama, D. W. 1968, Nature, 217, 511 
Sachs, R. K. & Wolfe, A. M. 1967, ApJ, 147, 73 

Smith, K. M. & Zaldarriaga, M. 201 1, Mon.Not.Roy.Astron.Soc, 417, 2 
Spergel, D. N. & Goldberg, D. M. 1999, Phys.Rev., D59, 103001 
Taylor A. & Watts, P 2001, Mon.Not.Roy.Astron.Soc, 328, 1027 
Verde, L. & Spergel, D. N. 2002, Phys.Rev., D65, 043007 
Yadav, A. P & Wandelt, B. D. 2010, Adv.Astron., 2010, 565248 



Appendix A: Tlie simulated L-ISW CMB bispectrum 
from the covariance method 

This appendix refers to the covariance method used to build the 
L-ISW simulated maps and described in Sec. 12.11 It is straight- 

z&„+Afl^- give the 



forward o check that the coefficients ^ 

JSW\3 



wanted bispectrum by calculating ((a^^^*^)-') — {{zcm + '^flf„,)"') 
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The lensing coefficients a^^^ can be expressed analytically, at first 
order in the lensing expansion, as 



«L + 



m+in'+m" y~>-mm'm" 



(A.l) 



{'{"m'm" 



C'({' + V)-{{{+\) + €"{€" + 1) 



where a^^ the primary and (p^^^ the harmonic coefficients of the 
lensing potential (p^. Since, according to the new variables defi- 
nition izcm, tcm) of eq.|6l tt„, 
p* 



A^^j, the fl^^j can be written as: 



ar.„ oc a 



4x10 



G 2x1 0~'" 



'x10~' 



'Cm ' 

Here, 



(A.2) 



for simplifying the notation, 

2 



ft 
SO 



that 



V ( ]\in+m' +m" /^-n\ 

lul't"m'm"\ ^> '^tCC' 

Aflf... - fca^*t(,„ at first order The explicit expression for 




1100 



(A3) 



{{a'^;/^^)^} takes the form: 

+ Z(,niAa'^„f + 3zlMn, + (A^nf + 2z(„,{Aaif} 
From this the only non-zero term is: 

{3z'cmA4j = 3</fflf;ff„,C[^(xf„,xp,„,(rf^)2+3;f,„3;f,„,,(l -(r[^)2) 
From this only survives: 

6{ftaZte,„Crxey^'^y(,,n' ^l-ir]^-). (A.5) 



1200 1300 
I 



1400 



1500 



By using the definition of in Eq. |2] and the approximation 
rj** << 1 for which 1 - (r^'^)^ - 1, since by construction: a^^^^ - 




recover the expected signal: 



1, {y ) - 1 and (xy) = we 



(A.6) 



Appendix B: L-ISW cross correlation coefficients 

The definition of the CMB lensing-ISW/RS cross correlation co- 
efficients is (TSperg el & Goldberg., 1999; .Verde & Spergel,2002"t 
iGiovi et al ] l20Q5h : 



Jo 



r(zis) - r{z) 



^P,p{k,z) 
oz 



dz. (B.l) 



k=- 



where, r{z) is the co-moving conformal distance and P^{k,z) 
is the gravitational potential power spectrum which accounts 
for both the linear and non-linear contributions. The non-linear 
regime RS contribution to the signa l is tiny, in agreement with 
(lLewisl2012tlJunk & Komatsul2012h Fl Considering both the lin- 
ear ISW and the Rees Sciama effect improves the /, 



NL 



vari- 



ance and signal to noise by few percent 2%) with respect to 



^ We found tha t the am plitude on the non-linea r effect estimated i n 
iMan gilli & Verde ( 2009') (as well as probably in iGiovi et al.l (120051) ') 
compared to Lewis ( 2012,): |junk & Komatsj j2012h was due to an inter- 
polation issue. Even a small numerical effect at the interpolation scale 
k = can propagate through the line-of-sight integration to a rele- 
vant effect on the non-linear transition scale of the C^**. Note however 
that the Cj^** coefficients are very sensitive to the cosmology parameters 
related to the late time evolution, t^ A, w, erg, and to the modeling of 
the non-linearities (e.g. lVerde & Sper gel (2002)), so extra care must be 
taken when comparing results from different authors. 



Fig. C.l. Noise simulations. The figure shows an example of 
one simulated anisotropic noise map realization (bottom panel) 
and its correspondent noise power spectrum in red on the upper 
panel. The black line corresponds to the power spectrum from 
the same L-ISW simulation. 



the linear only case calculation. In this work we consider both 
contributions for completeness. As a template, for both the sim- 
ulations and the estimator, we used the late ISW-lensing cross 
correlation coefficients of Eq. lB.ll This is a good approximation 
since this effect is the one which give s the main c ontribution. 
However, for a detailed description see (lLewisll2012b . 



Appendix C: The experimental setting 

Figs. IC.2] and IC. 3 1 summarize the experimental settings used in 
the simulations. These settings are inspired by a space-based 
experiment such as WMAP or Planck, with a variation of the 
noise with the ecliptic latitude. We consider a one channel CMB 
experiment with a Gaussian beam with a FWHM Qi, - 7', a 
galactic mask leaving ^ 80% of the sky and anisotropic un- 
correlated noise. In particular we obtain a galactic type mask 
from the IRASEI 100pm map smoothed at 5 angular degrees res- 
olution and with a threshold of 12MJy/sr. We consider a dipole 
type anisotropic noise covariance matrix, which accounts for the 
anisotropics owing to, e.g., the scanning strategy. An example 
of an anisotropic noise realization and the correspondent power 



http://www.cita.utoronto.ca/ mamd/IRIS/IrisOverview.html 
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9.0e-06 ^^^^^^K .^^^H 3.0e-05 ^^^^^^^^^^^^^H 13352 



Fig. C.2. Experimental setting: Beam window function and Fig. C.3. MoiA:. Galactic mask cut with /5t,, = 0.78 (upper panel) 

anisotropic noise map. The Gaussian beam with a 7' FWHM obtained from thresholding the smoothed lOO/zm IRAS map 

is shown in the upper panel, while the bottom panel shows the (bottom panel), 
dipole like anisotropic noise covariance matrix map. 

spectrum is given, respectively, in the bottom and in the upper 
panels of Fig. lCTI The noise start dominating from £ ^ 1300. 
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